Abstract 



We study the interactions of a Bragg-grating soliton with a locahzed defect 
which is a combined perturbation of the grating and refractive index. A fam- 
ily of exact analytical solutions for solitons trapped by the delta-like defect 
is found. Direct simulations demonstrate that, up to the numerical accuracy 
available, the trapped soliton is stable at a single value of its intrinsic pa- 
rameter 6. Depending on parameter values, simulations of collisions between 
moving solitons and the defect show that the soliton can get captured, pass 
through, or even bounce from the defect. If the defect is strong and the soli- 
ton is heavy enough, it may split into three fragments: trapped, transmitted, 
and reflected ones. 
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1 Introduction 



Optical fibers carrying a Bragg grating (BG) written on their cores have at- 
tracted a great deal of attention for various reasons [Q. Fibers equipped with 
homogeneous or chirped BG are a basis for the design of many elements used 
in optical telecommunications, such as filters p|, add-drop WDM multiplex- 
ers 1^, ^, dispersion compensators 0, and others. They also find important 
applications as sensors, since the optical properties of the fiber-based BG 
are strongly affected by local stress, see, e.g., a recent work P and reviews 
0, PJ. Besides that, BG-carrying fibers are interesting in their own right, as 
a unique physical medium supporting various optical excitations. In partic- 
ular, strong dispersion induced by BG opens a way to observe slow light in 
this medium, which is a topic of great current interest |]^, [TD[ . 

The BG-induced dispersion in combination with the Kerr nonlinearity of 
the fiber give rise to localized pulses, which are usually called gap solitons 
(or, more generically, BG solitons [|1T|). After being predicted analytically 
|1T2| , T3[| (see also the review |ll|]), BG solitons have been observed in ex- 
periments |T^ , They have attracted a lot of interest as their dynamical 
properties are very different, in many respects, from the properties of usual 
solitons in optical fibers (except for the regime in which the amplitude of 
the BG soliton is small, then it is virtually equivalent to a usual nonlinear- 
Schrodinger soliton ||16|). In particular, quest for very slow solitons, that are 
expected to exist in fiber BGs is a challenge [p!7 |. 



The subject of this work is interaction of a soliton with a local defect 
in the BG-carrying fiber. A local defect in this medium may trap light via 
four- wave mixing without the formation of a soliton per se [Q, but the 



interaction of solitons with defects is especially interesting |T^. First, it is a 
fundamental physical problem - in particular, because the defect may help 
to stop the soliton and thus to create a self-sustained pulse of standing light 
(leaving aside the problem of the fiber-loss compensation, which is a subject 
of another work |^). Second, narrow soliton pulses interacting with local 
inhomogeneities may increase the sensitivity when the grating is used as a 
sensor. Moreover, the trapping of soliton in a fiber BG may be used to design 
optical-memory devices. Further potential applications include optical delay 
lines and oscillators, see the concluding section of this paper. 

The objective of the present work is a detailed study of the interaction 
of BG solitons with a local point-like defect. In particular, it will be shown 
that a moving soliton can be captured by the defect. We will also study the 
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stability of trapped solitons, concluding that they may be both unstable and 
stable, depending on the values of the system's parameters. It is necessary 
to mention that the interaction between a BG soliton and local defect of 



the grating was considered in a very recent paper |T9[, which was chiefly 
focused on a detailed analysis of the case of finite-size defects. We, instead, 
consider point-like defects, which makes it possible to find an exact solution 
for trapped solitons, and an analytical potential of the interaction between 
a moving soliton and the defect. The study of the soliton-defect interaction 
reveals features that were not reported in Ref. - most notably, reversal 



of the sign of the interaction from attraction to repulsion in some cases. 
Another interesting property which we present here is the fact that a stable 
stationary trapped soliton matter of fact, an attractor: up to the 

accuracy of numerical simulations, it is found to exists at a single value of 
its energy, but initial solitons from a broad area in the system's parameter 
space relax to this stable state. 

The paper is organized as follows: in section 2, the model is formulated, 
and analytical results are presented in section 3. Section 4 contains numerical 
results for the stability of trapped solitons, which show the existence of the 
attractor. In section 5, simulations of collisions of moving solitons with the 
defect are presented, the most remarkable feature being the above-mentioned 
reversal of the sign of the interaction; this, in particular, gives rise to two 
different regions in the parameter space, where the moving soliton is captured 
or repelled by the defect. The paper is concluded by section 6. 



2 The model 



We approximate a local disruption (or enhancement) of the resonant Bragg 
grating due to the presence of a defect by a (5-function perturbation added 
to the standard nonlinear BG model for the right- and left-traveling waves 
u{x,t) and v{x,t). Taking into regard that the local perturbation may also 
affect the fiber's refractive index, cf. Ref. [M, we arrive at the following 



model, written in the usual notation 12, 13|: 



— h — h f + (-jMp + I^Hm = k5{x)v + X5{x)u, 

■9V .dv / 1 I |2 I l2\ e-/ \ A e-/ \ 

I— I— — Vu + [-\v\ + \u\ )v = Koixm + Aoixw, 

ot ox 2 



(1) 
(2) 
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where the ordinary ratio between the self-phase-modulation (SPM) and cross- 
phase-modulation (XPM) coefficients, 1 : 2, is adopted, the local defect is 
placed at the point a; = 0, and the real constants k > and A measure, 
respectively, the local suppression of the grating (k < corresponds to a 
locally enhanced Bragg reflection) and local perturbation of the refractive 
index. 

It should be mentioned that the region with a locally changed refractive 
index may also give rise to partial reflection, independently from that induced 
by BG. The corresponding perturbation is described by the cross-couphng 
linear terms proportional to the ^-function, which are similar to those (pro- 
portional to k) that are already present in Eqs. (|ID and (Q); however, the 
coefficients in front of these additional terms may be complex, due to the 
phase shift of the reflected wave. Thus, in the most general case, the terms 
n6{x)v and n6{x)u with real k in Eqs. (|I]) and (|^) can be replaced, respec- 
tively, by (ki + iK,2) S{x)v and (/ti — 1^2) 6{x)u, with two real parameters, ni 
and K,2- However, in this work we confine the consideration to the case of 
real k. As for the parameter A, it must be strictly real in the conservative 
model (its imaginary part accounts for a local gain, if it is present in the 
model [gg). 

If K is positive and small, the corresponding perturbation terms in Eqs. 
([1|) and (^ imply that the Bragg reflection is disrupted in a region whose 
size is much smaller than a characteristic length necessary for the complete 
reflection, which is, typically, ~ 1 mm [^, ^ , and which also determines 
a characteristic intrinsic spatial scale of the soliton. If, on the other hand, 
K ~ 1, this implies that the size of the disruption region is comparable to the 
reflection length. In the latter case, an exact form of the spatially distributed 
perturbation should be, strictly speaking, taken into regard (see Ref. 0]). 
Nevertheless, the use of the 5-function approximation may be justified in 
this case too, as the fundamental features of the interaction of a soliton 
with the inhomogeneity may be weakly sensitive to the exact shape of the 
inhomogeneity. The latter property is well known to be true in many other 
applications of the soliton perturbation theory [^. In accordance with this, 
in the numerical part of the work we will assume that k takes values between 
and 1. On the other hand, in the case k < 0, \k\ may take any value, as in 
this case the (5-function terms in Eqs. (0) and (H) have the literal meaning 
of a local perturbation in the form of the enhanced Bragg reflection. 

It is relevant to mention that this difference in the interpretation of 
positive and negative values of k, is very similar to the well-known differ- 
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ence between the positive and negative values of a perturbation parameter 
that determines the local suppression or enhancement (corresponding, re- 
spectively, to the so-called micro-resistor or micro-short) of the tunneling 
supercurrent in the standard model of a long Josephson junction, which is 
based on the sine-Gordon equation with a 5-functional perturbation |2T 



Indeed, a micro-resistor may only make the existing critical density of the 
supercurrent smaller, while a micro-short may arbitrarily increase the local 
value of the critical density. As concerns the refractive-index perturbation 
parameter A in Eqs. (||) and (0), it may take arbitrary values, as it literally 
accounts for a local perturbation. 

Notice that Eqs. (1) and (D can be derived from the Hamiltonian, 



H 



+00 



.du 



. dv du* 



ox ox ox 



dv* 
dx 



|4 I I |4 

u\ + \v\ 



I |2 I |2 

\u\ \v\ 



{u*v + UV*) 



In this expression, the asterisk stands for the complex conjugation, the overall 
sign of H is chosen so that the SPM and XPM terms correspond to self- 
focusing, and the perturbation part of the Hamiltonian, which accounts for 
the interaction of the soliton with the local defect, is 



vat 



k{u*v + uv*) + a (I 



|2 I 

u\ + 



\x=0 



(3) 



int ' 



3 Analytical results 

We start by seeking for stationary solutions to Eqs. (|1|) and (0) in the form 
u = U{x) exp {—iut) , V = V{x) exp {—iut) , (4) 
which leads to equations 



ljU + iU' + V + {(t\U\'^ + \V\^)U = n6{x)V + X6{x)U, (5) 
ujV ~iV' + U + {a\V\'^ + \U\'^)V = k6{x)U + X6{x)V, (6) 

the prime standing for d/dx. A family of exact solutions to the unperturbed 
equations (|^) and (||), with k = A = 0, is well known [0, |13| (these solutions 
are frequently called gap solitons): 



U{x) = J -(sin 



I sech ( (x 



- ■ sin - -9 



V 



-U*, 



(7) 
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where ^ is the coordinate of the sohton's center, and ^ is a real parameter 
taking values 

0<e <TT, (8) 
which also determines the frequency in Eqs. (H), 

uj = cos 9. (9) 

3.1 An exact solution for the pinned soliton 

An exact solution for a stationary soliton trapped by the local defect can be 
found in the form [cf. Eq. (|^] 



U = J -(sin 6*) sech 



1 

[x + a sgnx) ■ sin — -9 



, V = -U\ (10) 



where an expression for the parameter a is obtained by integrating Eqs. (^ 
and (P) in an infinitesimal vicinity of the point x = 0: 

/^T^tanh f - AV2) 
tanh (a sin ^) = W ^ — — — ^. (11) 

^ ' \ K + \ tan {9/2) ^ ' 

An example of this exact solution, with k = 0.4 and 9 = O.Svr, is shown in 
Fig. 1. 

In Eq. ([TT|), we assume that k"^ > (the opposite case will be considered 



separately below), and a/k^ — is realized in such a way that it has the same 
sign as k itself. In particular, in the case A = 0, Eq. ( pJ]) takes the form 

, , tanh(/s:/2) 
The intensity of the field at the central point of the trapped soliton is 



<Jsmh (asm 6'J + cos^ 

Note that |f/(a;)| is continuous in the solution ( p!OD across the point x = 0, 
while the phase of U{x) has a discontinuity at this point. 

If K < 0, Eqs. (|Tl|) and (|T2|) yield a < 0, hence, as it follows from Eq. 



2 sin^ 9 



(10), in this case |?7(x)| has two symmetric local maxima at a; = =Fa, and a 
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local minimum at the central point x = 0. On the contrary to this, if k > 0, 
we have a > 0, which implies that the solution has a single maximum at 
a; = 0. These features can be readily understood by dint of the interaction 
Hamiltonian (^. Indeed, together with the relation V = —U*, see Eq. (^, 
the expression yields the following value of the interaction Hamiltonian 
for the trapped soliton: 

(^intLapped = -2('^-A)|f/(x = 0)^ (14) 

With regard to the constraint /t^ > adopted above, we conclude that this 
contribution increases the net Hamiltonian if n is negative, and decreases it 
in the opposite case. Therefore, following the general principle stating that 
the ground state of a system corresponds to the minimum of H, the system 
tries to make \U {x = 0)f smaller if k < 0, and larger if k > 0. 

The result for the case < A^ can be obtained from Eq. (11) by analyt- 
ical continuation [or directly from Eqs. (H) and 



/X^^tan f VA2 - ^2/2) 
t^K(a=in9)^-^_^-^^, (15) 



where sgn (^v^A^ — = sgnA. In the particular case k = 0, Eq. ([15D 
simplifies to 

, / tan (A/2) , , 

cf. Eq. (p^ . Equations (|T5|) and (|16D yield a > 0, i.e., the single maximum 
of |f/(a;)| at X = 0, if A is negative, and a < 0, i.e., a local minimum at 
a; = and two maxima at a; = =Fa, if A is positive. Taking into regard 
the present constraint k"^ < A^, it is easy to check that these conclusions 
are again in accordance with the principle of minimization of the interaction 
Hamiltonian (|I1|). 

In the special case k = — A, both Eqs. ([TT|) and (|1^) amount to 

K A 

tanh (a sin 6*) = ,^ , , = ,^ , , . (17) 

^ ^ tan(^/2) tan(^/2) ^ ^ 

Lastly, in the case k = A, both expressions ([TTl) and ( |15|) yield a = 0, which 
is easy to understand: due to the relation V = —U*, the unperturbed soliton 
d^) remains an exact solution to Eqs. (^ and @ if k = X. 
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Obviously, the expressions (|n|) , (|T^ , and ( ]T3| ) through (|T^ make sense if 
they yield |tanh (asin^)| < 1, otherwise the solution for the trapped soliton 
does not exist. The latter restriction means that the exact solution exists not 
in the whole interval (H) in which the unperturbed solution was found, 
but rather in a narrower band. 



(18) 



where, in both cases /t^ > A^, 



2 tan 



' K 



tanh (Vfi;2 _ ^2/2) 



(19) 



This expression takes the form 6'niin = A in the case n = 0, and in the case 
A = —K it is 6'min = 2tan^^ (I'^D- 

The existence limit for the trapped solitons, as given by Eqs. (0) and 
([T9|), was also verified by numerical methods, solving Eqs. (^ and (P) with 
some approximations for 6{x), see details below. The numerical solutions (not 
shown here, as they are not especially interesting by themselves) demonstrate 
good agreement with the analytically predicted existence limit. 



3.2 A remark about the soliton stability 



Some (incomplete) conclusions about the stability of the trapped solitons 
can be made on the basis of the known Vakhitov-Kolokolov (VK) criterion 
22| , which states that a necessary condition of the stability of the soli- 
ton is dE/duj < 0, where oj is the soliton's frequency, see Eq. (^, and 
E = |w(x)|^ + |f (x)!^ dx is the energy (norm) of the solution. In the 

general case, the expression for the energy of the exact solution ([To|) is cum- 
bersome, but it takes a simple form in the case A = 0, on which the numerical 
simulations will be focused below: 



sin ^ ( sech k] 

2 ^ ^ 



sgn K 



(20) 



The calculation of the derivative dE/duj by means of Eqs. ([20| ) and 
d^) shows that all the trapped-soliton solutions satisfy the VK criterion. 
However, this result is not sufficient for claiming stability of the solutions. 
Indeed, in the unperturbed case (k = 0) the VK criterion per se predicts that 
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all the solitons (J^) are stable. However, studies of the stability, performed 
originally within the framework of the variational approximation ||23|, and 



then by means of accurate numerical methods |2^, had demonstrated that, 
in fact, the interval (P) in which the solitons exist is divided into two parts: 
they are stable if 

0<e<eer, ^cr~ 1.011 ■(7r/2), (21) 

and unstable if 6cr < 6 < vr. Notice that 6cr is just slightly larger than 7t/2, 
which is the value at which the gap sohton (0) has the largest amplitude and 
smallest width. This example clearly suggests that direct investigation of the 
solitons' stability is necessary in the present model too, which will be done 
by means of numerical methods in section 4. 

3.3 The potential of the sohton-defect interaction 

Besides the exact findings presented above, more general analytical results 
can be obtained in an approximate form, assuming that k and A are small 
parameters. In the lowest (adiabatic) approximation of the perturbation 



theory ||21| , one can find an effective potential of the interaction between the 
soliton, treated as a quasi-particle, and the local defect. To this end, the 
unperturbed wave form, based on Eqs. (|) and (0), is substituted into the 
expression ( |14|) for the interaction Hamiltonian, which yields 

" 3cosh(2esin^) + cos^ ^^^^ 

(recall is the coordinate of the soliton's center). In particular, in the case 
A = that will be considered below, the potential (^2]) clearly implies at- 
traction in the case k > 0, and repulsion in the opposite case. The latter 
circumstance strongly suggests that the trapped soliton cannot be stable in 
the case k < 0, > A^ fnor in the case A > 0, A^ > k^). 



4 Numerical simulations of the stability of 
the pinned soliton 

4.1 The approximation for the delta- function 

To simulate the stability of the trapped solitons and the interaction of a 
free soliton with the local defect, we have to adopt an approximation to 
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represent the 5-functions in Eqs. (^) and (^. In this work, we use a 

numerical scheme in which the coordinate x is represented by 501 grid points 
Xj, j = —250, —1, 0, +1, ... + 250. As an approximation to the 5-function, 
we have defined the following function on a set of 2N + 1 grid points in the 
central part of the integration domain, located symmetrically around zero: 

6 (a;„_(^^i)) ^ ^ (^TO^^) forn = 1, . . . , 2iV + 1, ^^3) 
^ ^0 otherwise. 

The normalization factor is defined so as to maintain the canonical normal- 
ization of the 5-function, 5{x)dx = J2j ^ i^j) = which yields 

-1 



A 




n-{N+l) ^ 
2N+1 



(24) 



Ax being the spacing of the grid (in fact, it was 0.04). 

In most cases presented below, we use N = 2 [then Eq. (p^ with Ax = 

0.04 yields A = ((l + v^) Ax) ~ 7. 726], which makes the (5-function quite 
narrow indeed. I3esides that, we also used another approximation for the 5- 
function, in which a simple rectangular regularization was used instead of 
(|23|). In some cases, sharp edges of the rectangular approximation produced 
problems with the convergence, but in those cases when results converged, 
they were virtually the same as generated by the approximation (pS]). 

It should be noted that, in the case kA > 1 [recall k is the perturbation 
parameter in Eqs. (|l]) and (H)], the local value of the coefficient k,6{x) may 
exceed 1, which is, strictly speaking, an unphysical situation, as it was dis- 
cussed above. However, we consider such a situation as modeling a physical 
one with a broader defect and smaller amplitude, when the integral strength 
of the perturbation, k S{x)dx, takes the same value, while the coefficient 
n6{x) does not exceed 1. Comparing results produced in simulations by the 
formally unphysical perturbation, and by its physical counterpart, we have 
verified that this assumption is, to a large extent, correct, with some minor 
differences observed only when the defect strength k is small. This issue will 
be further discussed below. 



4.2 The stationary soliton 

In simulations of the stability of the trapped solitons and the interaction of 
a free soliton with the local defect (the interaction is the subject of the next 
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section), we concentrated on the case A = 0, which is the most interesting 
one, as it focuses on the effects induced by the defect in the BG proper. We 
will mainly assume k > (although the opposite case will also be briefly 
considered), which corresponds to the most natural situation with the local 
perturbation suppressing the local Bragg reflectivity (rather than enhancing 
it). This is also justified because, as conjectured above and shown below, all 
the trapped solitons are unstable in the case n < 0. 

Our first objective is to analyze the stability of a stationary soliton 
trapped by the local defect. To this end, instead of using the exact so- 
lution given by Eqs. (p!OD and (|11]) through {^7^, we directly constructed 
trapped-soliton solutions to Eqs. (|^) and (H), with the regularized 5-function 
taken as per Eq. (^5]), by means of the Newton- Raphson method. An exam- 
ple of the thus found stationary soliton is shown in Fig. 1, along with the 
exact solution (0), for k = 0.4 and 9 = tt/2. After obtaining the stationary 
solutions in the numerical form, it was found that they are always close to 
the analytical solution for the ideal (5-function given by Eqs. (10) and (|TI1), 
so that the value of the parameter 9 corresponding to each numerically found 
solution could be easily identified. 



4.3 Stability simulations 

The stability of the stationary solutions was directly simulated by means 
of the split-step method applied to Eqs. (|I|) and (|^), employing the fast 
Fourier transform. First of all, in the case k < the numerical solution 
confirms that, as it follows from Eq. (|l^), |f^(a;)| has two symmetric local 
maxima at a; = =F a, and a local minimum at the central point a; = 0. 
Dynamical simulations confirm that all the trapped solitons with k < are 
unstable, as it was assumed above on the basis of the perturbation theory. 
The simulations further show that, as a result of the instability, two local 
maxima at x = =F a start to separate and eventually move away from the 
defect, so that the trapped soliton splits into two free separating ones. Since 
all the trapped solitons are unstable for k < 0, the investigation of this case 
was not carried out further. 

Typical results illustrating the stability of the trapped solitons in the case 
K > are displayed in Fig. 2, which shows the evolution of the initial pulse 
in the form of the trapped soliton (|l3), with its center placed at a; = 0, the 
value of the perturbation parameter being /t = 0.08 (note that this case is a 
"strictly physical" one, as it corresponds to nA < 1, see above). In the series 
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of panels displayed in Fig. 2, 9 is given values 6'in = OAtt, O.Svr, O.Tvr, O.Qvr. 
The left part of each panel shows a "side view" of the evolution of the soliton 
(in terms of values of \U\), the horizontal axis being time. The right parts of 
the panels (sometimes, these are insets) show a "top view" of the evolution, in 
terms of contour plots representing the solution. The evolution of solitons in 
the whole region k ^ 0.2 (note that it comprises both the "strictly physical" 
and formally unphysical situations) is found to be essentially the same as 
shown in Fig. 2 for k = 0.08. 

As it is seen from the plot in Fig. 2(b), and from the results for other 
values of k, only a single value of directly gives rise to a stable soliton (it is 
single up to the resolution provided by the accumulated data, i.e., we cannot 
rule out that it may be a narrow interval, rather than the single value). It 
is also found, more importantly, that in all the cases when the initial pulse 
relaxes into a stable soliton, the corresponding final value 6 stab of 9 is also 
a single one, being equal or very close to 7i/2, up to the accuracy of the 
numerical data. It is noted that 9 = -k /2 yields the maximum amplitude and 
minimum width of the soliton as per Eq. (p!0|), and it belongs to the stability 
region (^T]) of the unperturbed solitons, although being very close to the 
instability border. We also stress that the value 9 = 1^/2 always belongs, in 
the case A = 0, to the existence band (|T^) of the pinned-soliton solutions. 

The exact solution ( pUj ) for the trapped soliton with 9 = -n j2 takes (pro- 
vided that A = 0) a very simple form: as it follows from Eq. (0), in this 
case a = n/2. Note also that, if k is large, the latter result implies that the 
maximum value of the soliton's field is very small: according to Eq. (p!3D, 
\U{x = of = (4/3) (cosh/? — 1)~^. In fact, this means that, in the limit of 
large k, the trapped sohton becomes a quasi-linear defect-supported mode, 
cf. Ref. pf. 

The fact that there is a (presumably) single stable trapped soliton in the 
present model implies that the model gives rise to an effective attractor. 
Although attractors are usually considered as objects specific to dissipative 
systems, they may occur in conservative nonlinear-wave models due to the 
possibility of radiation losses. 

If is essentially smaller than ^stab, the soliton cannot relax to the stable 
state and may only decay into radiation, see Fig. 2(a). In fact, the expression 
(pop for the energy of the trapped soliton imposes an absolute limit on the 
values of 9^^^ of the initial unperturbed soliton (the one corresponding to 
K = 0) which can self-trap into the stable state. Indeed, the energy E^;^^ of 
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the initial unperturbed soliton is given by Eq. (|0D with k = 0, i.e., El^> = 
{8/3)6l^ , and, obviously, the energy Ef^^ of the established state can only be 
smaller than E[^. Taking for Ef[^ the expression (^OD with 6 = 6'stab ~ vr/2, 
we find i^fin = (8/3) sin~^ (sech/t). Thus, the limitation Efi-^ < e[^ takes a 
final form 

^J? > = sin-i (sech k) . (25) 

Obviously, {O^i^^ monotonically decreases with k, from 7r/2 at k = to 
0.22447r when k = 1 and further to zero as k — > oo. 

In order to investigate a possible dependence of 6'stab on the defect's 
strength k, simulations were run for a range of values of 9 slightly smaller 
than 7r/2, with k varying from O.lvr to O.Ovr. Similar to what is depicted 
in Fig. 2(a) and described above, solitons with 6'in < 7r/2 decay, but at a 
slower and slower rate, as is getting closer to 7r/2. For the simulation 
time t = SOOtt and k, = 0.1, the decrease of the amplitude of \u\ and |f | 
is larger than 1% if 9^^ < 0.497r. Of course, the decay stops as k ^ 0. 
However, the decay also becomes slower for larger values of n (which can be 
understood too: in the limit of k — >• +oo, the trapped soliton goes over into 
the above-mentioned linear defect mode, which is stable). For instance, at 
K = 0.9, the loss of the amplitude observed at t = SOOvr is less that 1% if 
0.477r < 6^in < 7r/2. So, there may be two possibilities: either 6'stab is getting 
slightly smaller than tt/2 when k is getting larger, or ^stab remains equal to 
71 / 2, but the solitons are less unstable, so the decay rate for 6'in < 6'stab is 
smaller for larger values of k. Because of the very slow decay or settling rate 
of |u| and \v\ for close to n/2, we cannot pinpoint a single value of 9 stab 
value; however, a very conservative estimate is that, at all the values of k 
considered, 9 stab is confined to the interval 0.477r < 9 stab < O.SStt. 

Figure 2(c) shows that if the initial value is somewhat larger than 
6'stab, the initial soliton evolves so that 9 relaxes to 6'stab, shedding off excess 
energy in the form of radiation waves. If is still larger for small values of 
K (actually, if 6'in ^ O.Svr at k = 0.08, see Fig. 3), the soliton does not relax 
into the stable state, but rather decays into radiation, sometimes generating 
a small-amplitude residual free soliton which is flung away from the point 
where the defect is located [see Fig. 2(d)]; spontaneous symmetry breaking, 
apparent in the latter case, is quite possible due to the instability. How- 
ever, for n ^ 0.2, the decay of the soliton at large values of 9 in is no longer 
observed. 

Figure 3 summarizes the results in the form of a stability region in the 
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plane (k, ^in), in which the initial soliton evolves into the stable one. It is 
quite natural that the stability region expands with the increase of k, as it 
is easier for a stronger defect to pull the field and thus capture the soliton. 
Besides that, Eq. ( pOD shows that the energy of the trapped soliton decreases 
with K, which also makes it easier to create a soliton from a given initial pulse. 
Lastly, one may notice that, if the value ^min limiting the existence region 



of the trapped solitons, see Eqs. (|18D and ([T9|) , exceeds 6'stab ~ 7r/2, the stable 
trapped soliton cannot exist. However, in the case under consideration, A = 
and K > 0, Eq. (|19D yields 6,nm = 2tan"^ (tanh(K/2)). Since tanh(K/2) < 1 
for all value of k, it is obvious that ^min is always smaller than 7r/2 (for 
example, 6*111111 = 0.2767r when k = 1), hence the condition Osmh > ^min never 
excludes 6'stab from the existence region. 



5 Trapping and reflection of a moving soliton 
by the localized defect 

The interaction of a moving soliton with the defect is a problem of obvious 
interest, both in its own right, and for such applications mentioned above as 
sensing by means of solitons and the soliton-based dynamical optical memory. 
To simulate the interactions, solutions for free moving solitons were first 
generated far from the localized defect in a numerical form by means of the 
Newton-Raphson method, as stationary solutions in the moving reference 



frame (these solutions are also available in the analytical form |]T2|, |T3| , which 
shows that the soliton's velocity c may take any value from —1 < c < +1). 
Then, full simulations of Eqs. (|lD and (H) were run by means of the split- 
step method combined with the fast Fourier transform, launching the moving 
soliton to collide it with the localized defect (it was again set A = and 
K > 0). 

First, the parameter plane (c, k) was explored, with c taking values 0.05, 0.075, 0.1, 0.2, 
and K ranging from 0.07 to 0.9. In all the simulations, the soliton was 
launched with 6 = 7r/2, which is equal (or very close) to the single value 
6'stab which, as it was demonstrated above, is selected by the localized defect 
as the single value at which the trapped soliton is stable. Figure 4 displays 
several different generic examples of the interaction of the moving soliton 
with the defect. In each panel, the lower part shows the evolution of the field 
|M(x,t)| by means of contour plots, and the upper part shows the waveforms 
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\u{x) \ and \v{x) \ at the end of the simulation, the location of the defect being 
marked by a short vertical line. 

If the defect is weak, the soliton passes through it, decreasing its velocity 
due to radiation loss. However, a stronger defect captures the soliton into an 
oscillatory state, see Fig. 4(a). The possibility of the capture of the soliton 
(through the radiation loss) is in agreement with the predicted attractive 
character of the interaction between the soliton and the defect, see Eq. (^21). 

Simulations demonstrate that, quite naturally, the minimum value of the 
defect's strength k, which is necessary for the capture of a given free soliton 
increases with its velocity. Figure 5 shows a border between the trapping (B) 
and transmission (A) regions in the plane (c, k), for fixed 6 = 7i/2. 

However, the soliton-defect interaction is more complicated. When the 
defect strength k is large, and when the soliton is energetic enough [either c or 
6 is large - region (C) in Fig. 5], its energy splits into three parts as a result of 
the collision: one part is trapped by the defect, and two others are scattered 
away, in both directions. The energy which is not trapped (either bouncing 
back or the passing through the defect) self-traps into a secondary soliton. 
Three examples of this mode of the interaction are shown in Figs. 4(c,d,e), 
for fixed values c = 0.5 and for 9 = O.Stt, and increasing values of the defect's 
strength, n = 0.3, 0.6, and 0.8. It can be seen that when n becomes larger, 
shares of the energy which are trapped and reflected back become larger. In 
regions A and B in Fig. 5, some energy is also scattered in the forward or 
backward directions; however, in these regions scattered energy is small, and 
it completely disperses into radiation, without giving rise to any secondary 
soliton. 

Another characteristic set of numerical data can be displayed for a fixed 
small value of the soliton's velocity, c = 0.075, while the parameter 6 of the 
moving soliton and the defect's strength k are varied. In this representation, 
the simulations demonstrate that the defect behaves in a qualitatively dif- 
ferent way for small and large values of 6. Namely, for 6 ^ 0.67r, the soliton 
interacts with the defect attractively, in accordance with the prediction of 
the perturbation theory, while for 6 ^ O.Gn the interaction becomes repulsive. 
The change of the sign of the interaction between the soliton and defect was 
inferred from the observation of outcomes of the collision. 

In the case 6 ^ 0.67r, the soliton gets captured by the defect, which clearly 
testifies that the interaction is attractive. The capture is accompanied by 
emission of a conspicuous jet of radiation. Passage of the soliton through 
the attractive defect is also possible (this happens, for instance, if the initial 
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velocity of the soliton increases). 

In a case when the defect strength is k = 0.5, and the sohton is heav- 
ier, having 9 = O.Tvr, the sohton bounces back from the defect, which is a 
characteristic feature of repulsion. The reversal of the interaction from the 
attraction to repulsion, which is observed for the solitons with 6 ^ O.Gvr, is 
beyond the applicability limits of the simple perturbation theory that gave 
rise to the effective interaction potential (|22|) . However, it agrees with the 
findings of the numerical stability analysis reported in the previous section, 
where it has been shown that solitons with large initial values of are un- 
stable if trapped by the defect, which may be realized as a result of repulsion. 

Figure 7 summarizes the results obtained for the interaction of the moving 
soliton with the local defect in the present case (fixed velocity). The border 
between attraction (A) and repulsion (B) regions is shown. To the left of 
the regions (A) and (B), there is a thin strip (C) where the defect is weak, 
K, ^0.1, and the soliton is heavy enough, 6 ^0A5. In this case, the soliton 
is neither captured nor reflected by the defect, but rather passes through 
it with significant radiation losses, see Fig. 6. By extrapolating the border 
between the regions (A) and (B), it can be inferred that the region (C) is 
made up of two sub-regions: an upper one, where the defect is effectively 
repulsive, and a lower one, where the defect is attractive. However, in either 
case the "heavy" soliton passes through the weak defect. 

Lastly, in the case when the defect is very strong, and the soliton is very 
heavy, the interaction may combine attractive and repulsive features. In 
this case, the energy of the incoming soliton is split into three parts: one 
is trapped at the defect, while the rest of the energy is scattered in both 
the forward and backward directions. The energies scattered away are large 
enough for self-trapping into secondary solitons. This happens in the region 
(D) of Fig. 7, and is similar to what happens in the region (C) of Fig. 5, 
where the defect was very strong, and the soliton was very fast. As the latter 
outcome of the collision is very similar to what was shown in Figs. 4(c,d,e), 
separate figures for this case are not included. 

As well as in Fig. 5, the region (D) in Fig. 7 is distinguished as a separate 
one to emphasize that significant amounts of energy are scattered away from 
the defect. In the other regions, some energy is scattered too, but its amount 
is much smaller. 

Concerning the results presented in this section, it is relevant to mention 
that moving unperturbed solitons also have their intrinsic stability limit. 



which, however, very weakly depends on the velocity |24]: as c increases 
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from to 1, the critical value Oct, above which the solitons are unstable, 
increases from 1.011 ■ (vr/2) [see Eq. (|2T|)] to 1.017- (vr/2). This implies that, 
strictly speaking, the incident soliton is unstable in all the cases shown in Fig. 
6; however, its instability is very weak, therefore it does not manifest itself 
in any way before the collision is completed (for the same reason, formally 
unstable solitons may seem fairly stable in an experiment). 

Finally, we return to the discussion of the validity of Eqs. (p3D and (p^ 
as an approximation for the defect when kA > 1, which corresponds to 
the formally unphysical situation. We stated above that such a situation is 
an appropriate model for a physical one, with a broader defect and smaller 
amplitude A, so that the product k6{x) does not exceed 1, provided that both 
the "physical" and "unphysical" perturbations are subject to the identical 
normalization, 5{x)dx = 1. 

As a possibility to verify this statement explicitly, we took the approx- 
imation (p3|) with = 19 [in the latter case, the normalization condition 
([2^ ) yields A = 1], and performed simulations of the collision of the soliton 
with this broad defect in the whole (k, 6) parameter plane, keeping the soli- 
ton's velocity constant, c = 0.075. The results were then compared to those 
displayed in Figure 7, which had been obtained for the narrow defect (recall 

= 2 in the case of Fig. 7). It was found that, to a large extent, the two 
sets of the results are similar indeed. The agreement is very good for larger 
values of k. Some disagreement occurs in the region (C), where k is small. In 
Fig. 7, the vertical part of the border of the region (C) is located at k ~ 0.1, 
with its tip at 6' ~ 0.45. For the broad defect with A^ = 19, it was found 
that this border moves out to k ^ 0.25, and the tip of the region moves 
to ~ 0.55, so that the region (C) is wider but shorter. This means that 
the conclusion whether the solitons can pass through the defect depends not 
solely on its integral strength k 6{x)dx, but also on its amplitude kA. 
However, since this affects a small part of the {k, 9) parameter plane, the 
above conjecture concerning the actual equivalence of the narrow (formally 
unphysical) and broad approximations for the (5-function is quite acceptable. 

6 Conclusion 

In this work, we have studied in detail the dynamics of the interaction of gap 
solitons with a localized defect, which is realized as a local suppression of the 
Bragg grating that supports the gap solitons. A family of exact solutions for 
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trapped solitons in the model with the 5 - functional defect, and an approx- 
imate expression for the effective potential of the soliton-defect interaction, 
were found in the analytical form. Direct simulations have demonstrated 
that this conservative model gives rise to an effective attractor: up to the 
accuracy provided by numerical data, the trapped soliton is stable at a single 
value of its intrinsic parameter ^ - in fact, at 9 ~ 7r/2, that corresponds to 
the soliton with the largest amplitude and smallest width. Depending on the 
strength k of the defect, initial solitons from a broad domain in the para- 
metric plane (k, 6) relax into the stable state; outside the attraction domain, 
they decay into radiation (sometimes, generating a residual free soliton). 

Simulations of the interaction between moving sohtons and the defect 
produce a number of different outcomes. If the soliton's velocity and 9 are not 
too large, the soliton is captured by the defect, in accordance with the fact the 
effective potential of the interaction between them is attractive. When the 
soliton is "heavy", with 9 ^ 0.6, the character of the interaction reverses from 
attraction to repulsion. In this case, the soliton bounces back. However, if the 
defect is weak, regardless of whether the interaction is attractive or repulsive, 
heavy solitons can pass through the defect. If it is strong, and the soliton is 
heavy or fast enough, the collision may split the soliton in three parts. Some 
energy is trapped by the defect to form a pinned soliton, while large amounts 
of energy scattered away may generate secondary solitons (transmitted and 
reflected ones) of signiflcant amplitude. 

These results suggest various applications. Besides the optical-memory 
design and sensing (defect detection), which were mentioned in the Introduc- 
tion, the delay of a soliton passing through the repulsive defect can be used 
to construct an optical delay line. Besides that, an optical resonator may 
be implemented by means of a soliton trapped in a cavity formed by two 
localized repulsive defects, with a locally enhanced Bragg reflection. In such 
setups, the local defects can be switched on and off in a controllable way by 
means of electrostriction. 
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Figure Captions 



Fig. 1. The real (even) and imaginary (odd) parts of the field U{x) for the 
trapped soliton, with k = 0.4 and 9 = O.Svr. The dashed curves depict the 
exact analytical solution (0) obtained for the ideal 5-function. The solid 
curves show the numerical solution, with the 6 - function approximated as 
per Eq. (p3| ) with N = 2. The two pairs of curves are nearly indiscernable. It 
was also checked that using an approximation closer to the ideal 6 - function 
brings the numerically found waveform still closer to the exact solution. 

Fig. 2. Typical results illustrating the stability and instability of the 
trapped solitons in the case A = and k = 0.08. The left part of each 
panel shows a "side view" of the evolution of \U\, starting with the initial 
stationary trapped soliton configuration. Right parts of the panels (or insets) 
show the evolution in terms of contour plots. T is the total simulation time, 
(a) 6*111 = 0.47r is smaller than ^stab- The sohton decays into radiation, (b) 
6*111 = O.Stt ^ 9 stab- The field directly self-traps into a stable sohton. (c) 
9 = O.Ttt > 6^stab- The field evolves into a stable soliton, shedding off excess 
energy in the form of radiation, (d) 9 = O.Qvr. This value is much larger than 
^stab, and the pulse decays into radiation, generating a small residual soliton 
thrown away from the defect. 

Fig. 3. The region in the plane (k, 6*111) between the two borders gives rise 
to the stable soliton trapped by the local defect with A = 0. In fact, in the 
region of relatively large values of k, the lower border becomes "fuzzy" : as it 
is explained in the text, no systematic deviation from 9 = n /2 was found, but 
a deviation within a margin of A9 = O.OSvr is not ruled out, as the instability 
of unstable solitons becomes very slow, impeding to monitor their evolution 
up to a definite result. Initial solitons taken outside the stability region decay 
(in the region above the upper stability border, a residual small-amplitude 
soliton can be found, being fiung away from the localized defect). 

Fig. 4. Collision between a moving soliton with the fixed value 9 = 7i/2 
and the defect. The lower and upper panels show, respectively, the evolution 
of the field \u{x,t)\, and the waveforms |it(a;)| (solid line) and \v{x)\ (dashed 
line) at the end of the simulation, (a) The defect with k = 0.2 captures 
the soliton with the velocity c = 0.2, which is accompanied by conspicuous 
emission of radiation, (b) In the case of a stronger defect, with k = 0.6, 
less energy passes through, while larger shares of the energy are trapped at 
the defect and refiected back, (c) The case of the strongest defect, with 
K = 0.8. In this case, most energy is refiected back, while the amounts of 
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energy trapped and transmitted through the defect arc smaller. 

Fig. 5. Regions in the parametric plane (c, k), with fixed 9 = 7i/2, in 
which the moving soliton, respectively, passes the defect, is captured by it, 
or is split into three parts. 

Fig. 6. The collision of the moving soliton with the defect. The parame- 
ters are c = 0.075, k — 0.1, 9 — O.Qtt. The upper and lower panels have the 
same meaning as in Fig. 4. In this case, the soliton overcomes the repulsive 
barrier and passes through the defect, with a considerable loss of energy. 

Fig. 7. A summary of results obtained for the interaction of the moving 
soliton and local defect, in the case when the soliton's velocity is kept constant 
at c = 0.075. In the region A, the defect is attractive for the solitons with 
small 9, and the solitons is captured by it. In the region B, the defect is 
repulsive for the solitons with larger 9; the soliton bounces back. In the region 
C, the defect is weak. The solitons passes through it, regardless of whether 
the interaction is attractive or repulsive. In the region D, the defect is strong, 
and, simultaneously, the soliton is heavy. The interaction combines repulsive 
and attractive features, and the soliton splits into three parts. One part 
is trapped, the remaining energy being scattered in forward and backward 
directions. 
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